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An optically pumped Cesium vapor magnetometer was used to 
measure the temporal fluctuations of the total geomagnetic 
field at Monterey, California. Power spectra were obtained 
from these recordings for the frequency range 0.1 - 10 Hz. 
Measurements were made during the month of April 1978 with the 
Kp index varying from 0+ - 8+. The power spectra obtained 
displayed a characteristic -20 dB/decade slope from 0.1 - 2.0 
Hz and were essentially flat between 2.0 - 10. Hz. Measurement 
times were: Local night (0000-0200), local morning (0800-1000) 
and local afternoon (1600-1800). Local morning exhibited the 
highest power density levels, exceeding the lowest power den- 
Sity levels (recorded during local night) by 10dB at .1 Hz and 
3dB at 10 Hz. A comparison of a magnetically disturbed day 
(Fredricksburg a index of 7) with a magnetically quiet day 
(Fredricksburg a index of 3) at the same local time (1730-1930) 
showed that at the low frequency end of the spectrum (.1 - .6 


Hz) there was 10dB more power on the disturbed day. 
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io eH NERODUCT ION 


It is often convenient to consider the geomagnetic field 
as made up of two major components: (1) The main field which 
can be considered to be approximately aumelle and generated in 
the interior of the earth, and whose observable temporal fluc- 
tuations are on the order of years or longer. For purposes of 
this study, we shall consider the main geomagnetic field to be 
constant in time. (2) A time-varying component superimposed 
on the main geomagnetic field whose origin is to be found in 
the ionosphere and the magnetosphere. This time-varying com- 
ponent of the geomagnetic field is of small amplitude compared 
femtne Main field (1 part in 10' or less) and involves frequency 
components up to the kHz range. 

The study of this time-varying eee of the earth's 
magnetic field is of interest in gaining a better understanding 
of basic geophysical processes as well as in improving the per- 
formance of existing or projected devices which are employed by 
the Navy. A number of operational Navy systems such as 
M(agnetic) A(nomaly) D(Cetection), magnetic mines, torpedoes and 
others, as well as some projected systems (such as ELF Communi- 
cation Systems), depend on the detection of small changes in 
the existing magnetic field in the ocean. Thus, from both a 
scientific and a Naval applications point of view, detailed 
experimental data on the geomagnetic field fluctuations seem 


highly desirable. 





As a first step in an experimental program to characterize 


and understand more fully the magnetic field 


LLUCHUa tlens ein 


and near the ocean, we have measured the power spectrum of 


mee iluctuations in Monterey, California using an optically 


pumped, total field magnetometer. -The measurements were 


carried out at a relatively quiet field site 
Village. A total of 52 hours of geomagnetic 
were made, at various times of day and under 


of geomagnetic activity. Frequency analysis 


near La Mesa 
field recordings 
varying conditions 


was performed on 


the data yielding the power spectrum of the fluctuations in 


the frequency range of 0.1 to 10.0 Hz. The general shape of 


the power spectra agrees with earlier measurements but several 


novel features were also observed. 
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II. BACKGROUND 


A. SOURCES OF GEOMAGNETIC FLUCTUATIONS 

Figure 1 shows an approximate power spectrum of the fluc- 
tuations measured on the earth's surface and having periods of 
a day or shorter. Also indicated on the graph are the approxi- 
mate amplitudes of the field fluctuations showing that even the 
largest of these fluctuations are less than 1 percent of the 
magnitude of the Main Field (B = 50,000nT). Nevertheless, they 
are of considerable interest both from the point of view of 
geophysics and from the point of view of Navy applications, 
and a world-wide network of magnetic observatories maintains 
continuous records of measurements of these varying field 
components. World Data Center A (operated by NOAA in Boulder, 
Colorado) issues several indices as well as forecasts of geo- 
Magnetic activity on a daily basis. 

The major diurnal variations have a period of 24 hours and 
are caused by current systems generated in the upper atmosphere 
as a result of solar heating. At mid-latitudes this Sq effect 
tends to be 25 - 50 nT, but in the equatorial regions, the 
effect of the equatorial electrojet can cause a typical vari- 
ation of about one hundred nT even during magnetically quiet 
conditions. There also exists a much smaller diurnal variation 
with a period of 12 hours, caused by lunar atmospheric tides. 
The magnitude of this variation is less than 0.1 of the daily 
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Figure 1. Power spectrum of geomagnetic disturbances observed on the 
Earth's surface. 


[Cladis, (1977) ] 
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Spherical harmonic analysis of these so-called Sq and L 
variations from a world network of observatories allows one to 
determine the equivalent electric current systems responsible 
for the observed magnetic variations. About two-thirds of the 
variations are found to be due to currents flowing external to 
the earth, and the remaining third can be attributed to currents 
induced in the earth's surface layers by the variable external 
field. The external currents are found to flow mostly at an 
altitude of about 100 km and are produced by convective move- 
ments of the conducting upper atmosphere across the earth's 
magnetic field lines. By analogy with the electric dynamo, 
this mechanism is often referred to as the atmospheric dynamo. 

Magnetic storms, which are severe, world-wide magnetic 
disturbances, have their main frequency components in the same 
region of the spectrum as diurnal variations. However, the 
rate of occurrence of magnetic storms seems to be primarily 
influenced by the solar sunspot cycle. 

Magnetograms obtained at low and mid-latitude stations 
indicate that many storms have a characteristic pattern of 
development as the storm progresses. The start of a typical 
storm, signalled by an abrupt world-wide increase in the hori- 
zontal component of the magnetic field, is known as a sudden 
commencement. This increase is typically 20 - 30 nT and has 
a rise time of a few minutes. The next 2 - 8 hours, during 
which the field remains at values higher than pre-storm, is 


known as the initial phase. This is followed by a 12 - 24 hour 


the! 





period during which the field decreases to values typically 
100 nT below its undisturbed magnitude, and this is known as 
the main phase. The final stage of the storm is characterized 
by a gradual return of -the field to its pre-storm value and 
lasts one to three days. 

Individual storm records will, of course, show deviations 
Peomethis “typical storm", particularly at higher latitudes, 
where rapid, large-scale fluctuations are common during storms. 
Even at mid-latitudes large amplitude fluctuations with periods 
of about one-half hour are occasionally observed during the 
main phase of a storm. 

Sudden changes in the dynamic pressure of the solar wind 
impinging on the outer boundary of the geomagnetic field are 
believed to be responsible for the itdetarslon of magnetic 
storms and the subsequent development represents the response 
of the magnetosphere and the charged particles within it to 
this initial impulse. 

The frequency region of greatest interest to us lies below 
3 Hz, commonly known as ULF. Geomagnetic field fluctuations 
in this frequency region are usually referred to as micro- 
pulsations, with periods varying from 0.2 - 600 seconds, and 
amplitudes from a fraction of an nT at the high frequencies to 
several tens of nT at the low frequency end of the interval. 
Micropulsations are divided into two general types: Continuous 
or pe and irregular or pi. The pe micropulsations display 


amplitude variations which are quasi-sinusoidal and last 
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several hours, whereas pi micropulsations have considerable 
irregularities both in frequency and amplitude. As indicated 
in Figure 1, the spectrum of micropulsations has been sub- 
divided into roughly 5 segments ranging from pe 1 at the high 
frequency end to pe 5 at the lowest frequencies. The occur- 
rence of pe's has a diurnal variation with a maximum near local 
noon. It is believed that pe's represent small hydromagnetic 
disturbances existing in the daytime magnetosphere, but the 
possibility that they are of lonospheric origin cannot be ruled 
out at this time. The pi's are trains of micropulsations con- 
Sisting of series of damped oscillations lasting from a few 
minutes to about one hour and having periods of roughly 40 - 100 
seconds. Their amplitude is largest in the auroral zone and 
decreases with decreasing latitude, and they appear to be pre- 
ferentially generated in the night-time magnetosphere. At high 
latitudes, there is evidence that pi's are hydromagnetic waves 
which have been greatly amplified by plasma instabilities 
created in the magnetosphere by high-energy electrons. 

In the ELF frequency range (3 - 3000 Hz), a number of 
emissions contribute to the magnetic fields, the most pronounced 
ones being the earth-ionosphere cavity resonances (Schumann 
resonances). These signals result from disturbances that are 
resonantly excited by lightning transients in the concentric 
Spherical cavity formed by the earth's surface and the lower 
region of the ionosphere. The power spectra of the ELF signals 
Often Show Maxima near 7.8, 14.1, 20.3, 26.4 and 32.5 Hz 


[Schumann and K6énig, 1954]. 
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Be REVIEW OF EARLIER WORK 

Power spectrum analysis of the geomagnetic field has been 
a popular tool to better understand geomagnetic field activity 
and to gain insight into possible source mechanisms for the 
fluctuations. For many naval applications one must also under- 
stand these geomagnetic phenomena, so that in attempting to 
detect man-made magnetic anomalies, one might try to eliminate 
this background signal. This allows for concentration on a 
smaller frequency range, which is of importance to the military 
observer. The next few pages are devoted to a short review of 
work done in power spectrum analysis of the geomagnetic field. 

In 1963, Santirocco and Parker reported [Santirocco and 
Parker 1961] their work on the polarization and power spectra 
of pe micropulsations. Recordings of the N-S, E-W, and verti- 
cal components of the geomagnetic field in the frequency band 
-005 - 5 Hz were made at a shore location in Bermuda. The field 
sensors utilized were ferromagnetic-cored loop antennas, polari- 


zed so that dB/dt vectors pointing North, East and down were 


+ Geomagnetic micropulsations are field variations with 
periods ranging from .2 seconds to approximately 10 minutes. 
A terminology recommended by the International Association of 
Geomagnetism and Aeronomy is pe for pulsations continuous and 
pi for pulsations irregular. The period ranges are: 


pe | .2 - 5 see Pier 1 - 40 sec 
pe 2 a eeOsee pare 40 -150 sec 
pe 3 10 - 45 sec 
pe 4 45 -150 sec 
peo 150 -600 sec 


For further details see Jacobs [1970], p. 16-20. 
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defined as positive. The sensitivity of the loops was 
70Luv-sec/nT] (where 1 nT = 1 Gamma or 10°° Gauss). With 350 
hours of continuous taping of the field components reported, 
power spectrum analysis of the tapes was carried out. The 
power spectra ofthe pe pulsations appeared as clumps of 
variable fine structure components superimposed on a residual 
background which falls monotonically with a slope of -6dB/ 
octave. Santirocco and Parker also observed that upon com- 
parison with the results of others that this background slope 
Memeonstant from .0002 Hz to nearly 1 Hz. 

Davidson reported in 1964 results on the average diurnal 
characteristics of geomagnetic power spectra in the period 
range 4.5 - 1000 seconds [Davidson 1964]. This study aimed 
at investigating the source and propagation mechanisms by 
means of the variability ag the power spectral density aa a 
function of 1) the spectral distribution of the source or 
driving function; 2) the coupling or impedance matching to the 
medium for various propagating modes, and 3) the propagation 
characteristics of the medium including the geometrical effects 
of the boundaries. 

The data was taken in the New Jersey coastal plain area 
(geomagnetic coordinates 39° 38.5'N , 74° 31.8'W and magnetic 
inclination 71°), using a Rubidium vapor magnetometer. Spectra 
were then obtained by digital processing, using the techniques 
of Blackman and Tukey [1968]. From these measurements, it was 


found that there are empirical expressions for the power density 


ay, 








Gow g 


ea) (in units on [ He) as a function of period as follows: 





for the 2.5 - 10.0 sec. period range the expression is 


e. 1.3 
er) = £226 X 10 * Tf (T in seconds). 


For the greater than 10.0 sec. period range the expression is 


GF) = 6.31 X Gites T*-© (T in seconds). 


It was also found that these two spectral components appeared 
to be independent of time. In addition, these period ranges 
contain spectral peaks superimposed on the constant slopes at 
9 sec. and 10 sec. periods. A second group of spectral peaks 
fall within the 17 - 70 sec. period range. 

Later another study made in the same area [Herron 1967] 
reported that in the 1700 - 2300 hour local time frame, the 
power spectra had the highest overall power density levels, and 
the strongest spectral peaks. It was also observed that the 
log log plot of Power vs. Period had a slope of -7.1 dB per 
Setave. 

Recently the efforts of A.C. Fraser-Smith and J.L. Buxton 
[1975], using a SQUID (Superconducting Quantum Interference 
Device) magnetometer:., enabled them to achieve higher sensitivity 
measurements of the geomagnetic activity in the frequency range 
-1 - 14 Hz. The measurements were taken at a mid-latitude loca- 
tion, in Stanford, California (geomagnetic latitude 43.5°N) on 
the N-S field component and spanned a 2-month interval between 


26 January and 26 March 1974. They found a minimum in the 
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spectrum in the 3 - 7 Hz frequency range while in the .1 - 3 
Hz range, there is a monotonic decrease with frequency. The 
variation is approximately f *7° for the .2 - 5 Hz (pe 1) range. 
For the frequency range 7 - 14 Hz, the magnetic activity is 
dominated by the first Schumann resonance. 

Finally, a report was made on the power spectrum observa- 
tions made from data collected in nine induction magnetometer 
observatories in widely separated locations about the earth 
[Wertz and Campbell 1976]. It was found that the slope of the 
log (Power Spectral Density) vs. log (Period) was about 3 at 
low latitudes and 4 in the Auroral Zone. There is a maximum 
of power spectral density at auroral latitudes and a minimum 
at the lower middle latitudes. For 10 second period measure- 
ments, the power spectral density varied from 2 X Ons EO 
2X10? (nT)*/Hz. For the 100 second period data, the power 
spectral density varied from 4 X 10°! to 2 X 10% (nT)?/Hz. The 
correlation between spectral density amplitudes and sunspot 
number was .9 for 12 second period data taken in College, 
Alaska. Earth surface conductivity, ionospheric variations 
and magnetospheric conditions were found to control the ampli- 
mide Of the signals. 

In addition to these published results, there is consider- 
able on-going work associated with various Navy projects, which 
involves some measurements of geomagnetic fluctuations in the 
frequency range considered in this report. For further details 


see Ref. 11, 13 and 14. 
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Pee ose oRUMENTAL SETUP 


A. GENERAL LAYOUT OF EXPERIMENT 

As mentioned earlier, there are essentially two components 
to the geomagnetic field. These are the main field and the 
time-varying components superimposed on the main field. This 
experiment dealt with the measurement of the time-varying field 
exclusively, treating the main field as the stable or d.c. com- 
ponent. 

The experiment was conducted in Monterey, California, in 
the community called La Mesa Village (Geomagnetic Latitude 
42.5°N). The sensors for magnetic field measurement were loca- 
ted in an adjacent wooded area subject to relatively few 
magnetic field disturbances. 

The measurement of the field was conducted using a Cs Vapor, 
optically pumped magnetometer with a pair of sensors positioned 
a distance of 55 meters (180 ft) apart on an East-West line. 

The main electronics and power supply equipment were housed in 
Building 318, 61 meters (200 ft) from either sensor. Electrical 
power to and signal from the sensors were transmitted via a 
coaxial cable 61 meters long. 

The geomagnetic field's time variations were measured at 
various times both local day and night, to get a sampling of the 
field as a function of time of day. The general level of geo- 


magnetic activity ranged from very quiet to disturbed during the 
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month of observations. Magnetic activity index Kp range was 


oF to st. 


B. DESCRIPTION OF MAGNETOMETER 

The sensing unit of the magnetometer is the Cesium Vapor 
cell. The operation of this sensing unit is based on the 
optical pumping and monitoring technique [Dehmelt 1957]. This 
technique has its basis in the Zeeman splitting of energy 
levels of the cesium valence electron in the presence of a 
Magnetic field. The energy levels are split into sublevels 
which have energy differences proportional to the total inten- 
Sity of the ambient magnetic field. 

To measure the energy differences between the sublevels we 
note that AE = hv and hence a measurement of v, the frequency 
of the emitted or absorbed photons determines AE which in turn 


1s proportional to the ambient magnetic field B since 


(g = Lande g factor, e = charge on the electron, h = Planck's 
constant, m = mass of the electron, B = the ambient magnetic 
field.) 

The measurement of v is accomplished by monitoring the 
optical transparency of the pumped Cs vapor cell while applying 
a RF field slowly swept in frequency around v. When the driving 
frequency corresponds to AE, the optically pumped electrons are 
redistributed among all of the sublevels leading to an increased 


absorption of the optical radiation. Thus at this particular 
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value of v, usually called the Larmor frequency Vy > we have 


GE = hy; - a B or Vy = rer B 

The constant of proportionality between the Larmor frequ- 
ency and the magnetic field intensity for cesium is 3.499 
Hz/nT. The cesium sensor actually detects the magnetic field 
strength by monitoring the transmission of resonance radiation 
through the optically pumped cesium vapor cell, and produces 
a Signal whose frequency is proportional to the total magnetic 
field intensity (with very fast response even to large field 
changes). Figure 2 illustrates the major components of the 
cesium sensor. A more detailed description of the optically 
pumped magnetometer is given in Bloom [1961]. 

The cesium sensor employed throughout this experiment was 
of the single cell design, which has large dead zones. This 
means that the sensor can only detect the magnetic field when 
the sensor axis is within the prescribed cone in relation to 
the total field vector to be measured. Figure 3 illustrates 
the active zones of the single cell sensor and shows that they 
include the greater part of one hemisphere. Optimum signal to 
noise ratio is achieved when the optical axis of the sensor is 
45° from the total field vector. Figure 4 shows the proper 
alignment of a signal cell sensor in the geomagnetic field, 


for maximum signal to noise considerations. 
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The sensor is effectively an oscillator which produces an 
rf signal (Larmor Frequency) directly proportional to the total 
field intensity. Thus the Larmor Frequency provides direct 
information on the field and can be used in conjunction with a 
Simple frequency counter. However, the more customary arrange- 
ment is shown in Figure 5, indicating that the Larmor Frequency 
signal 1S converted toa slowly varying voltage, whose ampli- 
tude is proportional to the frequency excursion of the Larmor 
Frequency, and hence the magnetic field at the sensor. 

Finally, the following procedures had to be followed for 
reliable magnetometer system operation: 

(1) The sensor was mounted in a relatively noise-free area, 
secluded from power lines, etc., and positioned on a firm non- 
magnetic support; 

(2) The sensor was securely fastened to the support to 
preclude its movement; 

(3) The cable end of the sensor was pointed in the magnetic 
North direction; 

(4) For optimum signal to noise considerations, the sensor 
axis was aligned at a 45° angle from the total field vector 
(Magnetic Latitude 42.5°N and magnetic inclination = 62° for 
Monterey, California); 

(5) The sensor was extended the full six feet from the 
sensor electronics box, and 

(6) The sensor and sensor electronics units were sealed in 
fiberglass housings to prevent water intrusion and provide shel- 


ter from the wind. 
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The specifications supplied for the Varian Model 4938 mag- 


netometer by Varian Associates are listed here: 


sensitivity T5 picotesla (.005 gamma) 
Information Rate Continuous and essentially instantaneous 
Range Continuous from 20,000-80,000 nT 


C. DATA COLLECTION SYSTEM 

The major components required to record the data on mag- 
Beere tape for analysis and display at a later time are: 

(1) Sensor; 

(2) Sensor Electronics; 

wey Readout; 

Geist terential Amplifier; 

So) Output Filter, and 

(6) Instrumentation Tape Recorder. 

Figure 6 shows the system components required to collect 
magnetic field data as an analog signal (i.e., voltage ampli- 
muae as a function of field). In this configuration, it is 
possible to tape record the analog signal for further analysis. 
Real time processing can be carried on simultaneously with tape 
recording. A comparison of the two spectra allows determination 
of the tape recorder noise. 

In this data acquisition configuration called the magneto- 
meter mode, the sensor generates the Larmor Frequency which is 
proportional to the ambient magnetic field. In the Readout, 


v 


the Larmor Frequency is mixed with a reference oscillator in 
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the mixer. The frequency of the reference oscillator cor- 
responds to the main field at the given location. In Monterey 
a reference oscillator frequency of 175,780 Hz corresponding to 
a 50.3 uT setting on the Readout was used. The mixer signal is 
a square wave having a frequency equal to the difference . 
between the Larmor Frequency and the reference oscillator 
frequency. This square wave is the input to the discriminator 
driver and discriminator, and thus the change in magnetic field 
is converted to a d.c. voltage whose amplitude is proportional 
to the magnetic field variation. In order to use the discrim- 
inator output, a differential amplifier is used to accept the 
discriminator's balanced output (balanced output has neither 
Side grounded). The output of the differential amplifier is 
then filtered (Bandlimited to .018-300 Hz) providing a suit- 
able signal for analog magnetic tape recording. Thus this 
system measures the variation in the total magnetic field inten- 
Smety . 

The signal was recorded on a Hewlett Packard 3960 instru- 
mentation tape recorder at a speed of 3 3/4 1.p.s. The flutter 
compensation option was employed to reduce the noise due to 
vibrations in and near the recorder. The signal sent to the 
recorder was a bandlimited (.018-300 Hz) analog output from the 
discriminator via the differential amplifier and the filter. 
The bandpass filter had slopes on low and high ends of 24dB/ 
octave. All data was recorded on Scotch Brand instrumentation 


tape, with each side holding 2 hours and 2 minutes of magnetic 
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field fluctuation data. All tapes were degaussed prior to 
using them in data collection. The data was recorded in the 
FM recording mode, which improves amplitude accuracy, and is 
less sensitive to tape imperfections. 

Some of the pertinent specifications for the tape recorder 


are listed below: 


Carrier center frequency at 3 3/4 i.p.s. Be oe chi 

Pass band =" oO Az 

S/N ratio without flutter compensation 4.8 

Distortion Less than 1.5% 

dee. drift 28 peak to peak 
output/degree 
Celsius 


D. SCHLUMBERGER MODEL 1510-03 SPECTRUM ANALYZER 

This spectrum er is a fully-digitalized high perform- 
ance instrument that provides real time spectrum component 
analysis of analog or digital input signals in 256 spectral 
bins across one of 10 selectable frequency bands. The result 
of each of the 256 independent analyses performed by the instru- 
ment is displayed on an internal C.R.O. as a dot indicating the 
power and frequency components of the analysis. The bandwidth 
of each line of resolution is dependent on the maximum frequency 
of the selected analysis range. The 1510 translates time domain 
Signals into the frequency domain by Fourier transform techni- 


ques. 


Syl 





There are five selectable input range attenuation switches 


werying from .1 Vea BO 10 oe full scale. In addition, the 


MS S 


1510 has power spectrum averaging selection controls to aver- 
age as many as 1024 individual spectra. 

The analyzer has an option of two window functions for 
application to the data in the time domain, as an amplitude 
function multiplied by the data samples. In effect, the entire 
analysis sequence is amplitude modulated by the window func- 
tion and the result transferred to the Digital Fourier Transform 
analysis memory. The window function is unity for the rectangu- 
lar, and for the Hamming function %-% cos 27 4 for each n 
of the 1 to 1024 data samples. The rectangular window has a 
smaller equivalent noise bandwidth and was thus employed in the 
data analyzing scheme. 

For the data analyzed in this experiment, the 1510 spectrum 
analyzer was used in one configuration only with the following 


settings kept constant for each sample: 


Input Attenuation eal! ae (OdB) 
Frequency Range le = 25.6 Hz 
Window Function Rectangular 
Mode Continuous 
Spectrum Axis _ X- log, f; Y-dB 
Stop after N depressed 
Averaging Interval SLY 
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fm CALIBRATION AND SYSTEM TESTS 

A number of system tests were performed prior to data 
collection. The first of these tests involved the frequency 
response of the magnetometer system. This test was designed to 
check the output of the discriminator of the readout by apply- 
ing aknown input to the sensor head. To accomplish this task 
the sensor and readout were set up aS illustrated in Figure 6. 
Measured field changes were produced and applied to the sensor 
by a sinusoidally stimulated 72 turn coil. A sinusoidally 
varying magnetic field component was generated by an a.ec. 
current of known amplitude and frequency flowing through the 
coil located as shown in Figure 7. 

The measured resistance of the coil was 1.1 ohms. The 
inductance of the coil was calculated from the equation for 


mranMetanece for a circular cross section coil. 


~ r2 8r : 
be N Ween (==) =>2) henries 
where N= number of turns = 72 
r = radius of coil = 33 X 10 *m 
a = radius of circular cross section = 1 X 10 *m 


Pee mecoledpi tty OL tree Space = 41 X 10772 


The inductance was calculated to be 7.7 mhenries. This set 
up a simple circuit for the stimulation process as shown in 
Figure 8. The impedance (X = RtwL) is a function of frequency 
(w) and the magnetic field generated by the coil is directly 


proportional to the magnitude of the current in the Gel wror 
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Figure 8 
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Magnetic Field Generation 
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the geometry described in Figure 7. Here the generated field 
lies along the earth's total field vector so that the coil's 
field either adds or subtracts from the total field vector. 


The magnetic field intensity along the coil axis is given by: 


fd 
H = =H es eed 


where N = number of turns 
I = eurrent in coil 


mes radlus Of coil 


im=mcistanece Of Coll to sensor 
Esliee = factor required for unit conversion from SLO lie a 
TG (one a meter 


jeer ieldk 


The coil-was "swung" with a sine wave of known frequency 
and amplitude to change the field a known amount. The field 
magnitude change was kept constant while varying the frequency. 
Figure 9 illustrates the results of this test. Clearly there 
appears to be a non-linearity of system response with frequency. 
The test was performed on two sensors with essentially identical 
results. It was noted that the frequency response was essen- 
tially flat from .1 to 2 Hz, for both sensors tested. The 
observed roll off from 2. to 10. Hz was approximately 1.1 dB 
per Hertz. The exact reason for this non-linearity is not 
known, but an appropriate correction for it was applied to the 
data. 

The next test involved determination of the noise level of 


the tape recorder that was used in data collection. This was 
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done in two ways. First a grounded input was applied to 
channel No. 1 of the tape recorder and the frequency spectrum 
was analyzed by the spectrum analyzer. Results are found in 
Figure 10. The analysis was run at tape recorder speed 3 3/4 
1.p.s. through a bandpass filter set for range .1 to 250 Hz. 
The spectra were averaged 256 times. Major noise components 
appeared at frequencies above the area of interest and were 
mostly power line (60 Hz) components or power line harmonics. 

Another test was used to check the reliability of the tape 
recorder to reproduce a given analog magnetic field signal. 
This test was designed to observe distortion of the signal by 
the tape recording process. Figure 11 gives the system compo- 
nents in block diagram form. A real time spectrum analysis 
-was run while the same signal was tape recorded. The tape 
recorded signal was then played back and also analyzed by the 
Spectrum analyzer. Graphs of the two spectra appear on Figure 
12. The tape recorded spectrum essentially duplicated the real 
time spectrum and the tape recorded curve is within £2 dB of 
the direct spectrum. The small spike at 3.3 Hz was later elimi- 
nated by using the flutter compensation feature of the tape 
recorder. 

It is not possible to put the cesium sensor into a field 
free region to measure sensor noise, since the sensor requires 
aeledst a 20.000 nT field to function. A readout noise measure- 
ment was possible however. This measurement was made prior to 


data taking, during data taking and again at its completion. 
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In order to obtain equipment noise associated with the read- 
out, filter and amplifiers, it is necessary to provide a 
constant frequency signal input as a replacement for the 
varying sensor signal. This constant frequency signal must 

be within the range of the total magnetic field present, so 
that the correct reference oscillator is selected. With the 
sensor disconnected, the synthesized "Larmor Frequency" was 
applied to the readout unit. The constant frequency synthe- 
sized "Larmor Signal" was beat against the reference oscillator 
in the mixer and the change in frequency is detected by the 
discriminator. Since the change in frequency is essentially 
zero, there should be no signal emanating from the discrimina- 
tor. This test will display any noise generated within the 
discriminator, mixer, filter and amplifiers. Figure 13 shows 
the results of an averaged noise spectrum taken over six periods 
within the span of one month. The noise spectrum is thus quite 
reproducible. The low frequency end of the noise spectrum 
exhibited the highest power content and dropped off by 20 dB 

ave 10 Hz. 

The power spectra obtained from the 1510 spectrum analyzer 
give only a relative dB scale on the y axis. It was necessary 
to calibrate the y axis to the power spectral density in y*/Hz 
or (nT)?/Hz. This was done following the procedure outlined 
in the EMR 1510 Spectrum Analyzer Operators Manual. The 
log, (nT)’ VS. log, (frequency) calibration is given here: 


lglvA 


. 
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Figure 13 
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Calibration 
Frequency Range iV COL or OE HZ 
Window Function (Rectangular) 0dB correction 
faput Signal 

(a) The calibration of the signal out of the discriminator 
and the differential amplifier was 4 mV per Hz of frequency 
excursion. - 

(b) Filter gain was accomplished by an amplifier within the 
filter and was a factor of 10.5 X. Thus the signal was 


10.5f4mV\_ 42mV 
Hz/ Hz 


(c) The cesium cell Larmor Frequency constant per nT of 


field intensity change is 3.499 Hz per nT, thus 


Wee, 2-p3dnz _ lu7mV¥ weference signal 


HZ gull gue 
STEP 1: Relate input reference signal to 100 mV, the 
OdB attenuation setting 


14 7mV =n )OOmVv = .68:) = ~-3.35dB: iSodk 
147mV 


Correction -3.35dB 
STEP 2: Read 3 correction factors from 1510 front panel 
PeEpeauency Pange .1 to 25.6 Hz - per Hertz cor- 
rection 
Correction +10dB 


b. Required attenuator setting to prevent overload 


Correction OdB 


ran 





c. Window function rectangular 
ENBW = 1 Eemmecciton —O0dp) 
bie. Add all Corrections in steps 1 to 3 for total 
correction 


= 300 dB 
mee. 00 SdB 


mon oods LOTAL CORRECTION 


This means that a signal corresponding to a 1lnT change in 


field strength at the 1510 input will require correction of 


Gees 
Az 


to be at the -6.65dB calibration mark. By relating the input 


meeede at 1510 output. Thus, on the C.R.O. is defined 
Signal to the 0dB level it is possible to calculate where Il1nT 
lies on the y axis of the dB calibrated C.R.O. The corrections 
are necessary to take into account what scales, window and 


frequency range are being used. 
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IV. EXPERIMENTAL RESULTS 


es INTRODUCTION 

In a series of preliminary tests and calibration runs, as 
eetianed in the previous section, the suitability of the mag- 
netometer and recording equipment to analyze geomagnetic 
fluctuations had been established. During the month of April, 
1978, a total of 53 hours of data was recorded, covering 
various times of day and a range of geomagnetic activity from 
very quiet to moderately disturbed. Each data run was 85 
minutes long from which 512 individual spectra were computed 
and then averaged. The width of each frequency "bin" was 0.1 
Hz or 256 bins between .1 and 25.6 Hz. The averaged estimate 
of the spectrum of a random input signal has a Chi-squared 
distribution, having two times the selected number of averaged 
spectra degrees of freedom. The statistical characteristics 
result in a 95 per cent confidence that the estimate will be 
within a ~.4dB accuracy for a 512 spectra average. 

The observations made of the magnetic field in Monterey 
were most concerned with the 0.1 to 10.0 Hz frequency range. 
Each of the spectra taken the spectral power density for this 
range was computed and displayed. Comparisons were made of 
the spectral power density at one specific time frame for a 
number of days. Comparisons were also made of spectra taken 
at different times of the day, and of the spectral power den- 


sities as affected by varying degress of field activity. 
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meee lyr TCAL SPECTRA 


Figures 14 through 18 each represent the average of 512 
Spectra taken from a recording of the earth's magnetic field 
at various time periods during the day. Each curve is a 
typical spectrum for daily field variations during the time 
frames indicated. 

Figures 19 through 21 illustrate the spectra of five 85 
minute runs for local night (0000-0200), local morning (0800- 
1000) and local afternoon (1600-1800), showing the day to day 
variations in the spectra. All times are given as Pacific 
Standard Time. 

Figures 22 through 24 are the averages of the five curves 
given in Figures 19, 20 and 21, except for local night, where 
only three of the five spectra were used. It was found that 
because one of the curves had significantly more power than 
the other four, and another of the curves had significantly 
less power, that it was decided not to include these two when 
computing the mean. Otherwise the high power spectrum would 
effectively be reproduced as the average spectrum. 

Figure 25 illustrates the power spectra of the 1730 to 
1930 period for an average day and for a magnetically disturbed 
day. The disturbed field has 10dB greater power level in the 
.1 to .6 Hz range. The spectra nearly duplicate each other 
from 3 Hz to 10 Hz. The Fredricksburg a index had a value of 7 
at the disturbed recording time and a value of 3 at the time 


the lower power curve was taken. These index readings were 
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provided by the National Geophysical and Solar-Terestial Data 


Center in Boulder, Colorado. 


me COMPOSITE SPECTRA 


Figure 26 is a composite of the important curves taken 
during the course of this experiment. It includes a five-day 
average for local morning and afternoon, and a three-day 
everage for local night. In addition, the equipment Bode 
spectrum is displayed. It is evident from this composite that 
the local night (0000-0200) spectrum approached the system 
noise Spectrum most nearly of the three periods observed. The 
local night spectrum had significantly less power than the 
other two periods only in the .1 to .4 Hz range. From .5 Hz 
to 10.0 Hz all three average spectra exceeded the system noise 
level by at least LOB. 

Local morning and afternoon had very similar spectra. At 
See terequencies the two spectra were separated by only 3dB 
in power. At all frequencies the morning spectrum had slightly 
more power than did the afternoon spectrum. Both spectra dis- 
played a definite -~20dB/decade slope in the 0.1 to 1.0 Hz range. 
Both curves also exhibited a drastic slope change from 1.0 Hz 
to 10 Hz. Here the slope changed to positive at 4 Hz and then 
back to negative at the 10 Hz point. 

The local night averaged curve had very low power at the 
low frequency end of the spectrum compared to the other two 
periods. This was true for the .1 to .4 Hz range. From .5 Hz 


to 10.0 Hz, the local night spectrum displayed the general 


60 





-2 F Figure 26 


Composite Spectrum 








Local Night (0000-0200) 
Local Morning (0800-1000) 
Local Afternoon (1600-1800) 


Noise of Equipment 


Power Dpectrol. Density \n (nt) /We 


edi 





trend exhibited by the local morning curve. In general, the 
local night curve had the lowest power density for all frequ- 
encies. 

The lowest signal to noise ratio occurred at the .3 Hz 
point where it had a value of 7.5dB. A signal to noise ratio 


of +10dB existed for all frequencies above .5 Hz. 
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Ve. CONCLUSIONS 


A. DATA SUMMARY AND RECOMMENDATIONS 

The characteristics of the spectra agree in general shape 
with those illustrated in the Trapped Radiation Handbook 
[Cladis, Davidson and Newkirk 1971], for the frequency range 
0.1 to 10.0 Hz. The power levels at these frequencies agree 
in general with those observed earlier at Stanford [Fraser- 
Smith and Buxton 1975]. 

The average spectra illustrate the varying level of magnetic 
fluctuations as a function of time of day. For the time periods 
analyzed, it was observed that local night (0000-0200) had the 
lowest power density for each frequency measured. The local 
morning (0800-1000) spectrum exhibited the greatest power 
density for each frequency measured. The .1 Hz point of the 
local morning spectrum had the greatest power density 
mee tOe-(nT)°/Hz). The 2 Hz point of the local night spectrum 
had the least power density recorded (4. X 107°(nT)*/Hz) thus 
spanning four orders of magnitude. 

It was found that the three average spectra converged as 
frequency increased toward 10 Hz, where only 2dB in power sepa- 
rated the night time and morning power density levels. 

During data collection, there were periods when the geo- 
magnetic activity level was too low to record. At these times 


only the equipment noise spectrum was observable. 
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The data taken during the experiment and presented above 
is deemed reliable, However, in order to extract meaningful 
geophysical conclusions, it is necessary to collect and 
Emalyze a far greater amount of data for different times of 
the year, different times of the day and for varying degrees 
of geomagnetic activity levels. It is therefore recommended 
that these measurements be continued and expanded, preferably 
with the addition of digitizing equipment for computerized 


data handling and analysis. 


B. EQUIPMENT IMPROVEMENTS RECOMMENDED 

The Varian 4938 magnetometer system provided reliable mag- 
netic field information in the 0.1 to 10.0 Hz frequency range. 
However, the frequency response of the system was not flat as 
Paeially anticipated, but rather rolled off by 8.5dB from 
2.0 to 10.0 Hz. The cause of this roll off is suspected to 
lie in the mixer-discriminator circuitry since both sensors 
exhibit essentially the same performance properties. 

The plotted power spectra did not show the same fine struc- 
ture as observed on the C.R.O. This is primarily due to the 
plotter-analyzer interface. Some fine structure elements were 
lost in transmitting the analyzer C.R.0. display to the plotter, 
due to a response mismatch between the analyzer and plotter. 

The single cell sensor was found to be reliable only when 
in a stable and stationary support. This type of sensor can 
not be allowed to have large angular motion, and must have 


stability in its orientation to the total field vector for 


optimum signal to noise characteristics. 
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The tape recording of the magnetometer data allowed for 
permanent record keeping of the data but contributed some 
noise to the system. The flutter compensation circuit cancel- 
led out much of the vibration noise and a good reproduction of 
the magnetometer discriminator output was thus achieved to 
Within * 2dB. A digital data collection scheme could greatly 
improve the absolute accuracy of the field measurements and 


eliminate the discriminator circuitry and its associated noise. 
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